phase portrait for non-linear system of DEs